S-Sirfc/- .SLui 


NATL AERONAUTICS AND SPACE ADM ; 


NA SA-C A-/Z¥/i>/ 




DYNAMIC ERRORS IN A TUNED FLEXURE-MOUNTED STRAPDOWN 
GYRO INTERIM SCIENTIFIC REPORT 

John E . Bortz 

Analytic Sciences Corporation 
Reading, Massachusetts 


September 1972 







^A-CR-12**nU DYWAHIC EChCrb Id a 
rtC TtlXUBE-nOOKED STEf.FDCMB GYBC 
tetio Scientific; Report. Jan. - 
72 (Analytic Sciences Corp.) t>* t o<- 
.25 


THt AWA1YTIC SCItMCH CO«fo <A TIOW 


TR-147-5 

dynamic errors in a tuner 
FLEXURE-MOUNTED STRAPDOWN GYRO 

Interim Scientific Report 
September 1972 


Prepared Under: 

Contract Number NAS 12-678 

for 

national aeronauts and space administration 

Marshall Space Flight Center 
Huntsville, Alabama 


Prepared by: 

John E . Bortz, Sr. 

Approved by: 

Arthur A. Sutherland, Jr 
Arthur Gelb 


the analytic sciences corporation 

6 Jacob Way 

Reading, Massachusetts 01867 

i 




f 


Mr. Ralph Kissel 
Technical Monitor 
NAS 12-678 

NASA Marshall Space Flight Center 
Huntsville, Alabama 35812 


j 

i 

j 

i 


a 


foreword 


MHe^f yt 8 L?« e,, | CeS Cor P° ration »»a B conducted a 

Flight c J^ 5 SUpp0rt 01 the Marshall Space 
Flight Center etrapdown technology program This 

1 mSZSSwF* d “ to « “• *™»t 

Mr. Jack H. Fagan of TASC performed the computer 
analysis work for this effort and Prof. Charles E 

the U f iversity of Massachusetts pro- 
vided extensive analytical support. Both contribu- 
tions are appreciatively acknowledged. 




ABSTRACT 


Motion induced errors in a tuned, flexure -mounted 
strapdown gyro are investigated. Analytic expres- 
sions are developed for errors induced by linear vi- 
brations, angular motion, and detuning. Sensor-level 
errors (gyro drift rate) and system-level errors (na- 
vigation errors) that are stimulated by an actual dy- 
namic motion environment are computed. 
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I 

1. INTRODUCTION ' 

I 

1.1 BACKGROUND | 

This report presents an analysis of the error torques acting I 

upon the rotor of a tuned, flexure -mounted gyro in a etrapdown applies- ' 

tion. This gyro was developed as a practical instrument in the early i 

19l »’s. Two of the first papers treating its development are by Howe 1 

(Ref. 1) and Savet (Ref. 2). A two-degree -of -freedom gyro, whose rotor i 

is supported by a gimbal configured as shown in Fig. 1.1-1, is an example 1 

of such a gyro. The rotor is attached to an internal gimbal by a flexure ; 

mount which ideally permits rotation only about one axis, as illustrated in 
the figure, and provides a restoring torque proportional to the angle of | 

rotation. The gimbal is connected to the spin motor shaft by another flex- 
ure mount, whose axis of flexure is nominally orthogonal to that of the | 

outer one and normal to the axis of the spin motor shaft, hi general, a 
tuned, flexure-mounted gyro can have any number of additional internal J 

gimbals functioning in parallel, with flexure mounts between them. 

In the configuration shown, the gyro can be arranged so there is J 

no average net torque on the rotor. This occurs because the dynamic mo- , 

tion of the gimbal, under a condition referred to as tuning (described anal- I 

ytically in Chapter 2), supplies an average inertia torque to the rotor that > 

exactly balances the average torque arising from the normal flexure of the * 

rotor supporting members (the flexure mounts). | 

Three types of^dynamic error; torque84n the tuned, flexure- v ^ 

mounted gyro are analyzed in this report: ! 


8 - 7 * 94 . 



Figure 1.1-1 Flexure -Mounted Gyro 


• Angular motion induced torques - 

• Specific force induced torques 

• Torques caused by rotor speed variations 

The errors generated by the rebalance loop itself, such as scale factor, 
errors and temperature sensitivities, are not analytically described al- 
though, in the simulation work performed in this^ study (refer to Chapter 5), 
scale factor assymetries have been included. 

In order to predict gyro drift under operational conditions, a 
dynamic environment has been chosen and the propagation of gyro errors 
torques into system level error in the presence of this dynamic environ- 
ment has been investigated. The description of the environment 
that is used in this study is obtained from a report (Ref. 3 ) which des- 
cribes the vibration spectra of the angular and linear vibrations and their 
cross correlations occurring in a CH-46C helicopter. In particular, 
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TASC’s System Error Prediction Program (SEP P)* l8 ueed to compute 

the system level errors generated by the Teledyne SDG-1 under the 
specified dynamic motion conditions. 

1.2 OVERVIEW OF BE PORT 

o* ^ The baSlC t0rqUC eqUatI ° nS for a flexure-mounted gyro are 
stated b, Chapter 2. The8 e stations are.then developed intone" 

form for the description of the dynamic errors in a practical tuned gyro. 

The dynamic response of a tuned, flexure -mounted gyro is des- 

Z J 1 , Z*” 3 ‘ The 6rr0r t0rqUeS lndUCed by case rotation are 

and t f , 3 ' 2 ' A fUnCtl0nal deSCrlpti0n of gyro is presented 

and trmisfer functions are derived. Consideration is given to the form and 

response of a general linear rebalance loop. Section 3.3 discusses the 

gular rate measurement process and the errors in that process. 

Chapter 4 is devoted to a mathematical derivation of the gyro 

errors due to case angular motion, case specific force, and rotor speed 
variations. 

The propagation of the dynamic errors into navigation errors 
as been studied for the case of a helicopter flight environment using the 
£*s em Error Prediction Program (SEPP). A description of the SEPP 
and the dynamic errors computed are presented in Chapter 5. 


* " “ 

This program is described in Ref, 4. 



1-3 


THE ANALYTIC SCIENCES CORPORATION 


EQUATIONS OF motion 

Equations are presented in this chapter describing the dynamics 

used^he 0 rT 8 “ Oti0n rClatiVe ‘° ^ CaSC - ^ C ° 0rdlnate 're- 
used in the development are introduced in Section 2. 1, the basic flexure- 

™ gyro equations are presented in Section 2. 2 , and in Section 2 3 

hese equations are put into practicai form for analysis of a tuned, ’ 

exure-mounted gyro for strapdown applications. 

2.1 COORDINATE FRAMES 

the fl TW ° al “ 0St COinCident coor dinate frames are used in analyzing 

■JTinrr b "°' °~ “ “• *•“ *. - «•* o. z 

T t “ or 

uie gyro case, about the C axis. Thus the £ and „ axes are in the 

otor plane of symmetry but do not spin with the rotor. This coordinate 

rame wi be called the rotor referenced frame or simply the rotor frame. 

axis and ^ 0,686 def * 10QD8 ’ * ^ 1,6 6eeD ^ Whe " the — shaft 
and rotor axis of symmetry are coincident, the C triad can be 

made to coincide with the x, y, z triad. m general, the case fixed and the 

rotor frames are not coincident and differ by the rotor deflection angles 

. * and e y , which define the rotation of the rotor frame relative to the case 

.™ e ' iS Shown m Flg - 2,1_1 - That is, 9 X represents a rotation of the 
rotor relative to the case about the x axis, etc. When 6 X and 6y are small, 
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R-73M 



V 


Figure 2.1-1 Gyro Coordinate Systems* 

as we shall assume they are, the two frames are related by the direction 
cosine matrix 



which transforms a case referenced vector x into the rotor coordinate 

— c 

frame according to the rule 

-r = C r -c 

• • • • . »: 

* ' 

A gyro with a single gimbal is shown for purposes of illustration. 

The development which follows is not restricted, however, to the 
specific configuration illustrated. 
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2.2 GYRO EQUATIONS: THE STATIC CASE 

A set of differential equations has been developed by Savet 
(Ref. 2) for the rotor deflection angles 0 x and 0 y for the static sit- 
uation in which there is no motion of the gyro case. These equa- 
tions are: 

( A+ !K +t V[ k -( a -!) w s] V (c+a > u a 6 ,* X9 y = -VOs2« s t-Vta2u. s t-M Jt -« x 

( 2 . 2 - 1 ) 

and 

( A+ !) § y +f ®y + [ k “( a '!) w s] 9 y ■ (C ‘ fa)a, s d x _X0 x = ^ cos 2^ ‘Sc® 1 " 2w s t ’ M y " € y 

( 2 . 2 - 2 ) 

where 

A,C(a, c) are the rotor (gimbal) moments of inertia about 
an axis in the plane normal to the axis of sym- 
metry and about the axis of symmetry respectively. 

0 X , 0 are the rotor deflection angles relative to 

^ the case about the case fixed x and y axes. 

f is the rotor deflection angle damping coeffi- 
cient. 

k is the torque constant for the gimbal and rotor 
flexure mounts. 

co g is the angular velocity of the rotor relative to the 
gyro case about the rotor axis of symmetry. 

X is a cross axis torque constant of the rotor shaft 
about the shaft axis of symmetry (for elaboration, 
see Ref. 2). 
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c x » *y are error torques. 

M x , M y a re torques applied externally to the rotor 
by the rotor torquing (rebalance) elec- 
tronics. 

The variables q x and q y are function of the rotor deflection angles and 
are given by 

^ = f S x +aw s *y -(*'!) "s®x (2.2-3) 

‘V = 2 8 y ' aa, 8 d x ' ( a ' I) "s ®y < 2 - 2-4) 

As can be seen from Eqs. (2.2-1) and (2.2-2), the flexure-mounted gyro 
exhibits the response characteristics of a coupled two^legree -of -freedom 
oscillator. The inertia coefficient A + a/2 is composed of a rotor con- 
tribution (A) and a gimbal contribution (a/2). A composite moment of 
inertia A can be defined as 

A = A+a/2 (2.2-5) 

The damping coefficient f is primarily attributable to uninten- 
tional phenomena such as viscous air damping. Damping may, of course 
be Intentionally introduced via the feedback torques M x and My, ’ 

The effective restraint torque coefficient K 

is of considerable interest in the flexure-mounted family of gyros. Neg- 
lecting, for the present, the drive shaft torque coefficient X, it can be 
seen from Eqs. (2.2-1) and (2.2-2) that if K = 0, then the gyro has the 
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response characteristics of a free gyro, i.e. , a gyro with no elastic 

restraint torques. 

For a given set of gimbal moments of inertia, a and c, and for 
a given flexure mount torque coefficient k, Eq. (2.2-6) shows that K is 
zero when 



Equation (2.2-7) is known as the tuning condition and a gyro for which 
Eq. (2.2-7) holds is known as a tuned gyro. As can be seen from the 
dependency of the tuning condition upon spin angular velocity, this 
is a dynamic tuning. It is made possible by the forced motion of the gim- 
bal, which exchanges energy with the rotor in an oscillatory manner. This 
energy exchange is phased so as to counteract, when averaged over one 

complete rotor spin cycle, the torque exerted on the rotor by the flexure 
mounts. 

2.3 A PRACTICAL TUNED GYRO 

An obvious simplification can be made in Eqs. (2.2-1) and (2. 2-2) 
for a tuned gyro, i.e., a gyro for which Eq. (2.2-7) holds. Also, for the pur- 
poses of analyzing the dominant dynamic error effects, it is reasonable 
to ignore friction effects in the flexure mounts and set f = 0. It is also 
reasonable to ignore as small the rotor shaft torques X6 and X6 With 
these simplifications, Eqs. (2.2-1) and (2. 2-2) reduce to 7 

AS x +C Vy = (2.3-1) 

A 0 y’ Cw s<§x = q y cos2 < o 8 t-q x sin2a )s t-M y -e y (2.3-2) 
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where 

C = C +a is the effective rotor moment of inertia 
about the axis of symmetry 

ax is the rotor angular velocity when the 
tuning condition is satisfied 

For a single gimbal configuration such as the one shown in 
Fig. 2.1-1, there can be an inherent stability problem. This arises 
if the rotor inertia moments satisfy the relationship C 2 = 2A . Then to 
the same degree of approximation, 

C = 2A (2.3-3) 

The homogeneous equation set related to Eqs. (2.3-1) and (2.3-2), is 
therefore 

X5 x + 2Aufe6 y = 0 (2.3-4) 

A0 y - 2A« s e x = 0 (2.3-5) 

The characteristic equation of this equation set is 

X 2 s 2 [s 2 +(2 ws) 2 J = 0 (2.3-6) 

where we have used Laplace operator notation. Equation (2.3-6) shows 
that the undamped natural frequency « n of. the system is 2n> s . As can be 
seen from Eqs. (2. 3-1) and (2.3-2), the forcing functions on the right 
haiid side are sinusoidal at a frequency 2o» g . This can cause unbounded • 
growth in 0 x and 
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The stability problem can be averted, as indicated by Craig . 
(Refs. 5, 6), by providing more than one gimbal and arranging the gim- 
bals in parallel. In the n gimbal configuration, there are n sets of os- 
cillatory forcing terms. . For n a 2 and for gimbal flexure mounts spaced 
2tr/n radians apart the forcing terms can be caused to sum to zero. Craig 
points out (Ref. 6) that "from a practical standpoint, at least three gim- 
bals are required for complete elimination" of the q x and q^ terms. 

It is> assumed that a multigimbal configuration is used to achieve 
stability and that the forcing terms that oscillate at 2u> g van is We 

therefore take for the basic equation set 


*8 x+ e* y 


- M x“x 


(2.3-7) 


AS -C6 

y * 


-M - e 

y y 


This set describes the dynamic performance of a tuned, flexure -mounted 
gyro whose case is not subjected to angular motion. Savet (Ref. 2) directs 
attention to the feet that, as indicated by the absence of spring restraint 
terms on the left hand sides of Eqs. (2.3-7) and (2.3-8), the tuned, 
flexure -mounted gyro exhibits the response characteristics of a free 
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3 * GYRO DYNAMIC RESPONSE 

The dynamic response of a tuned, flexure -mounted gyro is 
described in Chapter 2 for the situation in which the gyro case in not sub- 
jected to angular motion. In this chapter, we consider the general cir- 
cumstance of an arbitrary angular motion environment. 

3.1 ANGULAR MOTION DYNAMICS 

hi Appendix A, it is shown that the components of the rotor 
inertial reaction torque vector L are 

v' / v V 

L x = A ( A x +e x) +H (°y + e y ) + (c-A)n y n z 

+ e x [( C 'A)(°z * n y) + H« z ] 

+ e y [(C-A)nn y + (C-A) 0 2 ] (3.1-1) 

L y = A (‘"y + § y) ’ H (°x + 0 " (C ’ A) n x°z 
+e x [ (C ' A ) °x n y - (C - A) J 

" 6 y [(C ' A) (°z ' °x) ' H°z] (S- 1 - 2 ) 
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where 



°x’ ° z are the components of the angular velocity, 
of the gyro case with respect to inertial 
space 

® x > 0 are the deflection angles of the rotor rela- 
tive to the gyro case, as shown in Fig. 2. 1-1 

AfC are the rotor* moments of inertia about axes 
perpendicular and parallel respectively to 
the rotor axis of symmetry 

H = Cu; g magnitude of the rotor angular 

momentum 


As stated at the close of the preceding chapter, the tuned, 
flexure -mounted gyro Is essentially a free gyro. Therefore, the only 
external torques** applied to the gyro are those provided by the torque 
generator, M x and M y . Hence, to the degree of approximation for which 
Eqs. (3.1-1) and (3.1-2) are valid (see Appendix A), we can write 


M = -L 
x x 

(3.1-3) 

i 

ii 

sT 

(3.1-4) 


since the applied torques are ideally equal and opposite to the Inertia 
reaction torques. » 


The gimbal moments of inertia have been assumed to be zero. This as- 
sumption is based on the observation that although the gimbal inertias are 
of direct importance in tuning the gyro, they are of much lesser conse- 
quence in the dynamic motion error model, since AaA+a/2 and CaC-a. 

This assertion is relaxed in Section 4. 3 where the additional torques 
that exist when the gyro tuning condition of Eq. (2.2-7) is not met are 
considered. 
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Using Eqs. (3.1-3) and (3.1-4) in Eqs. (3.1-1) and (3.1-2) and 
rearranging slightly yields 

a 

A§ + Hd = -Aft -HO -M -e (3.1-5) 

x y x y x x ' ' 

AS -H$ - -Aft +HO -M -< (3.1-6) 

y x y x y y 



Eqs. (3.1-1) and (3.1-2) as 

<3t ang* Ai> * HC ' A) °A^ [ (C '^(“z ' °y) + H£ ^] + 9 y [ (C - A) n y °* + < C -*#,] 

(3.1-7) 

W AAy ' (C ' A) °A- 9 y[ (C - A) (°z -°x)- H °z] + e x[ (C - A) n x°y' (C ' A)ft z] 

(3.1-8) 

The subcript "ang" has been appended to the components and to 
Indicate that these are the error torques arising from gyro case angular 
motion. In Chapter 4 the error torques due to specific force and rotor 
speed variations will be discussed. Before proceeding with an examina- . 
ticn of the gyro error components, however, we first develop a functional 
description of the gyro and its rebalance loop. 

3.2 GYRO FUNCTIONAL DESCRIPTION 

rr r . . .*■ ■* .. 

Equations (3.1-5) and (3. i^) can be written using Laplace 
operator notation as 
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A * 2 Hs |N. -A* -H k r x k 


[-HS AsZj^J [ H -Asj|n y J [My] -[ <y J l '" 2 11 

where we have employed a vector -matrix format for convenience. Using 
the relation H ^2Atu g in Eq. (3.2-1) and inverting, we obtain 


'G(a) G(s) 1 (T Bb 

1^-; s \ 2^ 


-H O M € 

X X X 


0 G(s) G(s) 

. y . I 8 2w s. 


SJN wav 


where the function G(s) is 


G(s) = 


W : ) 


(3.2-2) 


(3.2-3) 


Equation (3.2-2) is diagrammed in Fig. 3.2-1. 

From the figure we see that if £ = M = 0, then, dropping the 
s subscript from oj g for convenience, 

•x ■ - a A^)(¥> %- - ® «, 


Using Eq. (3.2-3), this becomes 


8 = -ifl 

X S X 


(3.2-4) 


This is a particularly good approximation for the Teledyne SDG-1 gyro 
treated in Chapter 5. It results from the fact that, for that instrument, 
C ^ 2A. We will make the latter approximation in several places in this 
report without repeating the justification. 
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We conclude that in the presence of a steady input angular rate, a rebalance 
torque Is required that is proportional to that rate. 

Let us now consider a linear torquing scheme of the form 




(3.2-6) 


where Fj(b) is the direct feedback path transfer function and F (s) is the 
cross axis transfer functiai. Both are realizable linear networks. The 
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gyro and rebalance electrics then are configured as shown in Fig. 3 2-2 
Inserting Eq. (3.2-6) into Eq. (3.2-2) gives 

t f c (b)g(s) r d (s)G(8) f>)g<8) f„o)g<8)1 r 1 ^ r , r , 

~ TH3— 5 — T— + ®x , 9 iSl 

F c (sXi(s) F d (s)G(s) FJbXJ(b) FJs)C(b) „ « 

l s 5S— + -1_ jj^yj O y C£) Gfe) , y 

(3.2-7) 

which can be inverted to nht n <n 


. ^ j 8 + F « (8)G(8) *5= F c< a >°«» F c (b)C( 8) F d (8)G(8) 1 IV 

([ - F s («)G(s) +5 l F d (8)C(s) s + F d (s)C(8) t ^- F c (b)G( 6) Oy 


G(8, [lj + S- F c< s) ] -G(8)[s + L Fd ( a) ] 

G(a, [ 8 4 F d< 8 )] g < s >[3§4 f c ( s >] 


(3.2-8) 


Where Eq. (3.2-3) has been used for reduction of some of the expressions 
arising in the inversion and where 


A(s) * s 2 +G (s) -4 ~+4 F c (s) + 2sF d (s) + ^i (3.2-9) 
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Figure 3.2-2 Tuned, Flexure -Mounted Gyro 
and Rebalance Loops 


For the linear torquing scheme as represented by Eq. (3.2-6), 
when Eq. (3.2-8) is used, the torque M is found to be 


sF c (3)+ 511 G(s) ( F ^®) +I d( 8 >) 8F d < a)+G(9) ( F c (8, + F d (8) ) 

-flF d <s)-C(s) (fJ(»)+fJ(s)) sF c (s>+^G< 3 )(f^b) + F*<s))| 

F c ( ' h F d ( *>) ( f ! (,) 4 F d ( *>) V*>) 

sow (r c (.)-^j r d (.)) - *&.)(£ F c (.)*F d w).Sgl #.>) 


M 

1 


X 

■ 4^ 

I" 

M 

i 

I 

.y 

( 

\ 


(3.2-10) 


It is of interest to evaluate the steady-state torques generated by a con- 
stant angular velocity input vector O. . The final value theorem of Laplace 
transform theory states (Ref. 7 ) that if F(s) is the Laplace transform of 
f(t), then 

lim f(t) = lim sF(s) (3.2-11) 

t oo s-0 

To apply this theorem to Eq. (3.2-10) we need to make some assumptions 
about the transfer functions F d (s) and F c (s). .Assume that 


F d (0) = k d F c (0) = k c 

From Eq. (3.2-3), note that 


G(0) = | 


Using Eqs. (3.2-12) and (3.2-13) in Eq. (3.2-9), 


(3.2-12) 


(3.2-13) 


A(0) 


k 2 + k 2 
K d K c 


(3.2-14) 


The Laplace transforms of the constant input rate vectors 0 


and e are 


i 



e 

J 

1 

x ss 

1 

x ss 

| °< 8 > = a 

O 

1(8) = s 

c 

(O 

to 


y ss 


(3.2-15) 


Then using Eqs. (3.2-12) through (3.2-15) inEq. (3.2-10) and applying 
(3.2-11), we obtain 
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. yssj 

' o -h] rn 
h o o 


' 

X 8S 


€ 

X SB 

y 8 sj 


€ 

. %s. 


(3.2-16) 


Thus in the case of a constant angular velocity about the x axis, there is 
a steady-state torque applied about the y axis to balance the input moment. 
As is expected, the x axis balance torque is zero when there is no input 
rate about the y axis nor error torque about the x axis. 


3.3 ANGULAR RATE MEASUREMENTS 

In a strapdown gyro, the calibrated feedback torques are used 
as a measure of motion about the case input axes. For a linear rebalance 

loop, the measurement M is taken to be linearly related to the input an- 
angular rate vector. 


a xx a xy I [° x 


M Iv v K 


(3.3-1) 


where O is a vector of the estimated angular rates about the two sensitive 
axes and the a y of the geometry matrix have yet to be specified. B we 
assume that the error torques are zero and that a steady-state condition 
exists, then, using Eq. (3.2-1), we see that 


3-9 


THE ANALYTIC SCIENCES CORPORATION 


0 -H In 


N [ H °Jl°yJ 


(3.3-2) 


or, solving for ri x and O y , 


O = i M 
x H y 


y h M x 


(3.3-3) 


To investigate the errors in this measurement scheme we 
mus consider the general case where the error torques are not zero and 

Eq reaCted ' SUbStitUtlneE<1 - (3.3-2) into 

fcq. (3.2-1) yields 


As Hs 0 


-As -H 0 -hITo 1 


(3.3-4) 


[-Hs As*J[*J [ H -AsJ^J [ H oILl | fy 


If we define the rate measurement error vector 60 to be 


aA = A - o 


(3.3-5) 


thenEq. (3.3-4) becomes 


<n * . “ “IN f" °1K] M 


•» °J j [-hs [. „j[ s j [, y 


(3.3-6) 
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or, solving for 60, 



Equation (3.3-7) is the equation for the dynamic rate measurement errors 
in a two -degree -of -freedom strapdown gyro. 

Summary , — The operation of the flexure -mounted strapdown 
' gyro hiis been discussed in this chapter. The gyro dynamics equation was 
restructured to reflect the effects of gyro case angular velocity. 
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4. DYNAMIC ERROR MODEL 

A principal objective of this report is the development of a 
model for the dynamic errors of a tuned, flexure -mounted strapdownrgyro. . 
In this chapter, models are presented for dynamic errors due to angular 
motion, specific force, and rotor speed variations. 


4.1 ANGULAR MOTION ERROR TORQUES 

The dynamic response of the tuned, flexure -mounted gyro, when 
subjected to arbitrary angular motion, is treated in Chapter 3 from an 
instrument dynamics point of view. In this section the angular motion 
error torques, as presented in Eqs. (3.1-7) and (3.1-8), are discussed 
term by term. These equations are repeated here for convenience. 

\ng = (C " A) [( c -A)O x Oy +( e-A)Aj . 

< 4 -i -if. 

^ang = 

(4.1-2) 
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4.1,1 Anisoinertia 


The terms (C-A) ° y n z Eq. (4.1-1) and (C-A) in 
Eq. (4.1-2) are referred to as the anisoinertia error torques. They are 
critical terms because they can generate constant error torques (bias 
drift rates) when case angular rates are constant or when properly phased 
oscillatory angular motion is experienced. Note that the primary gyroscopic 
input for the x axis is HO = C u> O . The x-axis anisoinertia term is 

X s X 

(C-A) O H a: CO z O v /2. Therefore, the anisoinertia error torques may 
z y j 

be written 


e 


‘0 0 

x anisoinertia 

C 

y z 

€ 

2 

-0 0 

y an isoinertia 


X z 


However, such a balance is more difficult to achieve since it must be 
accomplished by variation of the rotor parameters alone. * For practical 
reasons, such a balance is not attained. 

In the event of oscillatory motion about the z axis (the spin 
axis), a phenomenon known as decoupling comes into play . This occurs 
when the rotor speed control loop stiffness tends to degrade. As oscilla- 
tion frequency about the spin axis begins to exceed the rotor speed control 
loop bandwidth, sinusoidal variations begin to occur in the speed of the 
rotor relative to the gyro case; the rotor spins at a constant rate relative 
to inertial space but does not follow the case oscillations about the spin 
axis. Consequently the inertia of the rotor about its axis of symmetry, 
C, is eliminated from the effective anisoinertia. For high frequency 
angular rates about the gyro spin axis Eq. (4.1-3) becomes 


We see, by comparing the terms in Eq. (4.1-3) with the input 

terms -Cw O and Cw O , that the ratio of the coefficient of the aniso- 
s y s x* 

inertia term to that of the input term is O z /2w s . For a constant non- 
negligible case angular velocity about the z axis, this could result in a 
significant error.* In a single degree -of-freedom gyro, an attempt is 
made to balance the moments of inertia in order to suppress the aniso- 
inertia term. This balancing is facilitated by the fact that in the single- 
degree -of-freedom gyro, the inertia moments C and A are moments for 
the entire suspended element which consists of the rotor and its rela- 
tively substantial supporting gimbal. Such a balance of inertias would 
also be desirable in a flexure -mounted gyro. In fact, as can be seen from 
Eqs. (4.1-1) and (4.1-2), this would eliminate many of the dynamic error 
torques due to angular motion as well as the anisoinertia error torque. 

♦ 

Quantitative estimates of this error are given in Chapter 5. 


€ x 


o n ' 

anisoinertia 

■ (f-c) 

y z 

c 


-n n 

y an iso inertia 


X Z 


(4.1-4) 


That is, for the case where C ^ 2A, the sign of the anisoinertia coefficient 
is reversed when O z occurs at frequencies above the rotor speed control 
bandwidth. 


To remove the (C-A) terms in Eqs. (4.1-1) and (4.1-2) would require 
a rotor shape approaching that of a sphere. 
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4.1.2 Cross Coupling 

The terms 0 X [(C-A)(n| - O^) + HO z ] in Eq. (4.1-1) and 
-6y[(C-A)(n2 _ q2j _ HO z ] in Eq. (4.1-2) are referred to as the cross - 
coupling error torques. Now 


H = Cw ; 
s 


C-A a C/2 


(4.1-5) 


and since a) s » Oj, i = x, y, z, the terms involving H predominate. We 
can write the dominant cross coupling terms as 


*cc 

f *cc 


HO 


HO_ 



re 


X 


e 

J 

yj 


(4.1-6) 


Using Eq. (3.2-8) and neglecting error torques we can relate the rotor 
deflection angles to the input rates in the form 


re i 


X 


e 

= - 

yj 

L 



fn 


X 


n 


L yj 


where 


Q^s) Q^s) 
-Q2< s ) Qj(s) 


_ ^F d (s)G(s) + ^.F c (s)G(s) 
Q l (8) Z{5) 

F c(s)GC)-^F d (.)0(.) 

Q 2 (s) 


(4.1-7) 


(4.1-8) 


(4.1-9) 


Combining Eq. (4.1-7) with Eq. (4.1-6), we obtain for the cross coupling 
torques 
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J 


fx CC 


O 

i 

- 


Q^s) Qgfs) 


“I 

r H o n ' 


X z 

i 

HO 0 

J 

y Z J 


(4.1-10) 


In order to estimate the significance of these terms, suppose 
that Q x , fly, and f> z are oscillatory at a frequency below the first break 
frequency in Qj^s) and Q 2 (s). In that case we can replace Qj(s) and Q 2 (s) 
by Q^(0) and respectively. With the aid of Eqs. (3.2-13) and (3.2-14), 
we see from Eqs. (4.1-8) and (4.1-9) that, ' f 


Hk. 


Qi(°) = 

1 v? 


H* 

<¥ 0) = T*- 

d c 


(4.1-11) 


(4.1-12) 


In a typical linear rebalance loop, at low frequencies the direct rebalance 
torque is much larger than the cross axis rebalance torque, i.e., kj » k c . 


Hence Eqs. (4.1-11) and (4.1-12) can be approximated by 

H A 


and 


Qi (°) " t ~ \ 


Q 2 (0) = 0 


(4.1-13) 


(4.1-14) 


Thus to a good approximation, the cross coupling torques of Eq. (4.1-10) 
are 
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-0 0 

*CC x z 

= k H 

Cy q -n n 

Xcc y z 


(4.1-15) 


These terms are potentially serious because of the rectification that 

occurs when 0 or 0 are sinusoids that are in phase with 0 . The 
x y z 

cross coupling error mechanism is also excited by constant input rate 

components 0 and O or O and O . 

x z y z 


4.1.3 Rotor -to-Case Misalignment 

The terms 9 y [(C-AJOyOx+fC-A)^] in Eq. (4.1-1) and 
e x [(C -A) ttyOx " ( c ’ A )^z ] in Eq. (4.1-2) arise from rotor -to-case dynamic 
misalignments. In a single-degree-of -freedom gyro (Ref. 8), such terms 
arise because of suspension system compliance or because of misalign- 
ments generated during the gyro manufacturing process. In a two-degree- 
of -freedom gyro, the rotor misalignments 6 y and 0 X arise as a consequence 
of the normal rotor deflections that occur when there is case rotation. 

In order to estimate the relative contribution of the rotor -to-case 
misaUgnment torques £ mig , consider the first term of each torque com- 
ponent. 


< C - A|I VMK 


(c-A)n n 

' y X 


(4.1-16) 


For comparison purposes, assume C a 2A and rewrite the coefficient ex- 
pression as follows: 
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(C-A)O O = s 0 n 
y x 2o> y * 


By inserting Eqs. (4.1-13) and (4.1-14) into Eq. (4.1-7), at frequencies 
below the first break frequency, the rotor deflection angles can be 
approximated by 


Equations (4.1-16) through (4.1-18) can be combined to obtain 


0 O 

xy H 

q /n ' 


y x \2 a; 


By comparing Eq. (4.1-19) with Eq. (4. 1-15), it is seen that the misalign- 
ment torque is similar in form to the cross coupling torque except that 
the misaUgnment torque components include the multipUcative factors 
^y/2u> s and O x /2co s for the x and y components respectively. Since 
fty» 0 x <<: 2tdg, it is reasoned that the misalignment torques are corres- 
pondingly less significant than the cross coupling torques. This conclu- 
sion is borne out by the results of the error propagation study described 
in Chapter 5. 
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4.1.4 Angular Acceleration Sensitivity 


The terms -A^ x in Eq. (3.1-5) and -Afty in Eq. (3.1-6), while 
not included in the definition of the angular motion error torques, can 
generate sizeable vibration induced errors at the system level. Equation 
(3.3-7) identifies the same terms as "dynamic” rate measurement errors 
(the second term on the right-hand side). In a vibration environment os- 
cillatory angular rates imply oscillatory angular accelerations, and, from 
Eq. (3.3-7), oscillatory errors in the angular rate measured by the gyro. 
These can combine with the true rates to rectify in the strapdown attitude 

♦♦ v - ' * * m t * 

computations. 

Suppose an angular oscillation takes place about the x-axis of a 
two-degree -of -freedom strapdown gyro. The frequency is assumed to be 
well within the gyro bandwidth. The indicated outputs of the gyro will be 
ft x = + 6 ^x 311(1 ^y = 6 ^y* From Eq. (3. 3-7) one component of 6fty will 

beAOx/H. In particular, if 


O = W sin vt 
x 


60 = - A ^ cos vt + . . . 

y h 


The angular error in coordinate frame orientation caused by 60^ will be 
e y = J«O y dt 


sin vt + . 
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' "‘A™ 2 [ww^^] 


(4.2-2) 


m is the rotor mass 

6 Z is the z component of the displacement 
of the rotor center of mass from the 
center of support* 

f x ; f y- f -z are the com Ponents of specific force 

K are the elements of the rotor compliance 
J matrix.** 


4.2.1 Mass Unbalance 


The mass unbalance error torques are 


= -m6 f 


(4.2-3) 


f y = m V„ 

y unbal * x 


(4.2-4) 


in . vibratory «w™»»t, thus, te™, avar.ga to Wb. tbo spe . 

” " “* a b~fB.«!=ponant, the average ntaae 

unbalance torques are also nm^Sro. 


mount axes. (For a multi-gimbal gyro with the \SnSSL » ,_8h ** 
parallel, an neaare «oant^. U fLte r « “at aS‘, ^ 

.■a™“he^^*a“ 1 S£^te™ o?7 “• samt 

y 3 sm 6 le degree -of-freedom gyro (see Ref. 8 ). 
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4.2.2 Compliance 

The compliance torques are given by 

2 f„ o a 


[VA-VA*^ - «*£ * fee- K„) Va] 


% =m fi Kf f -K ] 
y comp L *yzy zx 


f x y +K xz f f ' K zx x + ( K xx' K zz) f zx] 


The compliance terms (also referred to as the -g squared- terms) pro- 
duce non-zero average torques in a zero average vibration field because 
of the rectification that occurs. For example, assume that all cross 
compliances K lj( i / j, are zero and that 


* = - f sin ut 

j z 


Where v is the linear vibration frequency. Then 


(4.2-7) 


l 'compj 

This torque has an average value 


K - K 1 
zz yy 

(mf) ^ k sin^ vt 

xx ' zz 


€ 

x 


'K 

- K 

comp 

= | (mf) 2 

zz 

yy 

c 


K 

- K 

comp 


XX 

zz 


(4.2-8) 
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2 2 2 
co = co. + 2 6 coco, + 6co 
s t t 


(4.3-5) 


B 6 o) is small compared to toj, as will be assumed here, the last term 
can be ignored. Combining Eqs. (4. 3-1) through (4. 3-5) gives the de- 
tuning error torque components 


= -2 6 co 


detune 


'fF\ ), 


(4.3-6) 


€ = -2 6 to 

y detune 


(4.3-7) 


Since rotor speed can be controlled very precise 1$ the principal 
causes of detuning are likely to be errors in measuring the rotor inertias, 
a and c, and the flexure spring constant, k, or a shift in the value of the 
latter. Consequently, detuning is typically a quasi-static phenomenon. 
Using this point of view, Eqs. (4.3-6) and (4.3-7) are rewritten as 


(4.3-8) 


where k J . is the quasi-static detuning factor 
detune 


“detune = 2 W k l a -fJ 


(4.3-9) 


At constant or low frequency input rates, Eq. (4.1-18) can be used in 
Eq. (4.3-8) to obtain 
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For a tight rebalance loop and accurate tuning and will both 

be small. Consequently the detuning torques will likewise be small. It 
is noted that when the flexure mounted gyro is not used as a strapdown 
instrument, and 9^ may not necessarily be small and the detuning error 
torque may have correspondingly greater significance. 


4.4 SUMMARY 

Error torques excited by angular motion, specific force, and 
rotor speed variations have been discussed in this chapter. The most 
significant term for each torque mechanism is tabulated in Table 4.4-1. 


TABLE 4.4-1 

ERROR TORQUE SUMMARY 


Kamo 

x Component 

y Component 

Approximation 
Valid for 

Full Expression 
Given In 

AnisoincrtU 
(Low Frequency) 

£ o o 
5 yx 

- £ n 0 

2 x * 


Low frequencies 

Section 4.1.1 

An Iso Inertia 
(High Frequency) 

- y n n 

2 y x 

£ n n 

2 xx 


High frequencies 

Section 4.1.1 

Cross Coupling 

-k 0 0 

q * 

-k h n n 
q y * 


Low frequencies 

- Section 4.1.2 

Angular Acceleration 
Sensitivity 

An 

-V 

AO 

X 

~ H 


All frequencies 

Section 4.1.4 

Misalignment 

n 

-k h n n — 
q x y 2« fi 

-k H 0 0 

q y * 

n 

X 

2w s 

Low frequencies 

Section 4.1.3 

Mass Unbalance 


mi f 
X X 


All frequencies 

Section 4.2.1 

Compliance 

m 2 (K«-K„) I y f i 

m *<K»-K«>y. 


All frequencies 

Section 4.2.2 

Detuning 

k k n 

detune q x 

k k n 

detune q y 


Low frequencies 

Section 4.1.1 
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5. SYSTEM ERROR GENERATION 

In this chapter, the characteristics of the Teledyne SDG-1 
flexure -mounted gyro (Ref. 9), used as a strapdown sensor, are inserted 
into the TASC System Error Prediction Program (SEPP) to study motion - 
induced errors in that instrument. The SEPP is briefly described in 
Section 5.1 . In Section 5. 2 we discuss the results of the SEPP study, 
namely the dynamic errors of the gyro when subjected to a realistic dy- 
namic motion environment. A brief comparison of the performance of 
the SDG-1 gyro with that of another gyro in the same environment is also 
presented. 

5.1 SYSTEM ERROR PREDICTION PROGRAM 

The System Error Prediction Program is used to compute sys- 
tem level attitude drift rates and acceleration measurement errors result- 
ing from the interaction of strapdown IMU gyro error mechanisms with 
angular and linear vibrations. The information flow in the SEPP program 
is illustrated in Fig. 5.1-1. The motion environment can be described in 
detail by density spectra which relate the magnitude of the in-phase and 
quadrature components of all possible pairs of angular and linear vibra- 
tional motions in three dimensional space. Thirty six unique spectral 
density functions exist, (hi the usual situation where such a mass of data 
is not available, a reduced number of spectra can be used. ) The gyro and 
accelerometer error parameters (moments of inertia, rms scale factor 
errors, etc. ) and the sensor rebalance loop transfer functions are included 
in the error computations. The program does not treat nonlinear rebalance 
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5. SYSTEM ERROR GENERATION 
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x detune 
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y detune 
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]h this chapter, the characteristics of the Teledyne SDG-1 
flexure -mounted gyro (Ref. 9), used as a strapdown sensor, are inserted 
into the TASC System Error Prediction Program (SEPP) to study motion - 
induced errors in that instrument. The SEPP is briefly described in 
Section 5.1. In Section 5. 2 we discuss the results of the SEPP study, 
namely the dynamic errors of the gyro when subjected to a realistic dy- 
namic motion environment. A brief comparison of the performance of 
the SDG-1 gyro with that of another gyro in the same environment is also 
presented. 

5.1 SYSTEM ERROR PREDICTION PROGRAM 


For a tight rebalance loop and accurate tuning k^ and kjctune both 

be small. Consequently the detuning torques will likewise be small. It 

is noted that when the flexure mounted gyro is not used as a strapdown 

instrument, 0 and 6 may not necessarily be small and the detuning error 
x y 

torque may have correspondingly greater significance. ... - *. 

4.4 SUMMARY 

Error torques excited by angular motion, specific force, and 
rotor speed variations have been discussed in this chapter. The most 
significant term for each torque mechanism is tabulated in Table 4.4-1. 


The System Error Prediction Program is used to compute sys- 
tem level attitude drift rates and acceleration measurement errors result- 
ing from the interaction of strapdown IMU gyro error mechanisms with 
angular and linear vibrations. The information flow in the SEPP program 
is illustrated in Fig. 5.1-1. The motion environment can be described in 
detail by density spectra which relate the magnitude of the in-phase and 
quadrature components of all possible pairs of angular and linear vibra- 
tional motions in three dimensional space. Thirty six unique spectral 
density functions exist. (In the usual situation where such a mass of data 
is not available, a reduced number of spectra can be used. ) The gyro and 
1 accelerometer^error.paxameters (moments of inertia, rms scale factor 
errors, etc. ) and the sensor rebalance loop transfer functions are included 
in the error computations. The program does not treat nonlinear rebalance 


TABLE 4.4-1 

ERROR TORQUE SUMMARY 


Name 

x Component 

y Component 

Approximation 
Valid for 

Full Expression 
Given in 

An iso inertia 
(Low Frequency) 

sV. 

- ? n n 


Low frequencies 

Section 4.1. X 

An iso inertia 
(High Frequency) 

-§«A 

f 0 n 

2 x z 


High frequencies 

Section 4.1.1 

Cross Coupling 

-k n o 
q x z 

-kHOO 
q y z 


Low frequencies 

Section 4.1.2 

Angular Acceleration 
Sensitivity 

AO 

"H 1 

AO 

X 

IT 


All frequencies 

Section 4.1.4 

Misalignment 

n 

-k H n n 5 -*- 
q X y 2 w 8 

-k H no 
q y x 

0 

X 

2wg 

Low frequencies 

Section 4.1.3 

Mass Unbalance 

- m6 z f y 

m 6 f 
z X 


Alt frequencies 

Section 4.2.1 

Compliance 

“ W*yy) y. 



Ail frequencies . 

Section 4.2.2 

Detuning 

k detune k q °x 

k detune k q *y 


Low frequencies 

Section 4.1.1 
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loops in detail, but several pulse torqued loops can be approximated by 
equivalent linear continuous mechanizations which do meet the descrip- 
tion required by the program (see Ref. 4). The SEPP program gives a 
detailed source -by-source breakdown of sensor-level and system-level 
attitude drift rate and errors in measuring specific force. It allows the 
major causes of vibration -induced errors to be quickly identified and 
provides an indication of how well specific systems can do in various dy- 
motion environments. 


5.2 STRAPDOWN SYSTEM ERRORS 

The vibration power spectra used to represent the dynamic 
environment for the SEPP run is the helicopter environment described in 
Ref. 3. The Teledyne SDG-1 gyro rebalance loop characteristics are 
described in Ref. 9. Additional gyro parameters used in the runs a^e 
listed in Table 5.2-1. The accelerometers were assumed to be perfect 
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TABLE 5.2-1 

SDG-1 PERFORMANCE PARAMETERS 



m this analysis.** Coordinate frames were employed as shown in Fig. 5. 2-1. 
The SEPP runs were performed using two SDG-1 gyros with "Gyro 1" in- 
strumenting the x and y axes and "Gyro 2" instrumenting the z axis of the 
helicopter coordinate frame. Both sensor-level errors (gyro drift rate)and 
system-level errors (drift in the computed coordinate frame) and errors in 
the processed specific force * arising from the sensitivity of the SDG-1 gyro 
to the CH-46C helicopter dynamic environment^ were studied. 


Obtained in a telephone conversation on December 27, 1972 with 
0.02 S;g/hr/g2 edyne SyStemS C ° mpan y- Corresponds to 

tribute le hf ti0 L err f rS i° wWch ^Perfections in accelerometer con- 

!ta?assumZ e th t0 ** orally small in similar analyses 

that assume the same environment (see Ref. 4). * 

talS'SS force Itot toe been trace- 

,ra “ “* com Pensated /or angular 

n^Ii^Il e . UCOpter environment was chosen since it is one of the few re- 
presentative environments for which detailed angular and linear vibra 

Als °’ * has to ana[yze 

gyros, thereby providing a basis for comparison. 
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d ~ HELICOPTER DATA FRAME 

1 ~ GYRO 1 FRAME 

2 - GYRO 2 FRAME 

Figure 5.2-1 Coordinate Frames Employed in SEPP 

At the sensor level, four error mechanisms were found to give 
errors in excess of 0.001 degree per hour. These mechanisms are: 

• An isoinertia 

• Cross coupling 

• Angular acceleration sensitivity 

• Rebalance loop asymmetry 

Discussion of each of the mechanisms, with the exception of the rebalance 
loop assymetry, are presented in Chapter 4. 

The cross coupling error is analogous to that found in a single - 
4 • r degree -of -freedom, gyro, with one ‘ exception :’ lt can be generated by analog 
rates about either gyro input axis. Because of the normal deflection of the 
rotor spin axis in the operation of the gyro, the actual input axes are not 
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aligned with the reference input axes, which are fixed in the gyro case. 
Hence the gyro can also sense a component of angular velocity about the 
nominal spin axis. This sensitivity, which is treated analytically in . 
Section 4.1.2 is referred to as the cross coupling error mechanism. 

Angular acceleration sensitivity, which is discussed in Section 
4.1.4, is similar to output axis sensitivity in single -degree -of-freedom 
gyros. K is based on the fact that angular acceleration of the gyro case 
will cause output angles to grow. This growth is indistinguishable from 
that caused by the angular rates the gyro seeks to^measure. In particu- 
lar, an angular acceleration about one input axis has the same effect as 
an angular rate about the other input axis. Furthermore the error pro- 
duced is often of just the right phase relative to true angular motion to 
produce a rectification in the system attitude calculations. 

Rebalance loop assymetry refers to the difference between the 
rebalance scale factor when torquing in a positive sense and the scale 
factor when torquing in a negative sense. 

The gyro drift rate generated by the interaction of the helicopter 
dynamic environment with each of the error mechanisms listed above is 
shown graphically in Fig. 5.2-2. , As can be seen from the difference in — 
error components, the propagation characteristics are sensitive to gyro 
orientation. 

In addition to the constant gyro drift rates induced by the dy- 
namic motion environment, there is also induced a class of errors referred 
to as " system -level” errors. There are system-level drifts known as un- 
detected coning and pseudo-coning and system-level acceleration errors 
known as undetected sculling and pseudo sculling. Basically, these errors 
arise from the failure of the gyros (and accelerometers also in the general 
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case where perfect accelerometers have not been postulated) to detect 
high frequency oscillations or because oscillatory sensor errors are 
rectified in the system-level calculations. (For a treatment of the sys- 
tem level error phenomena, the reader is directed to Ref. 4. ) 

The system-level errors induced by the helicopter environ- 
ment are shown graphically in Figs. 5.2-3 and 5.2-4. The only error 
mechanism that contributes in any meaningful way to either pseudo-coning 
or pseudo-sculling error is the cross axis sensitivity of the instrument to 
angular acceleration (analogous to the output axis sensitivity in a single- 
degree-of -freedom gyro). Undetected coning and sculling errors arise 
because the dynamic motion spectral bandwidth for the helicopter exceeds 
the measurement bandwidth of the gyros and hence there are high fre- 
quency angular oscillations not recognized by the gyro. 



Figure 5.2-3 Vibration Induced System Drift Rates 



Figure 5.2-4 Vibration Induced Pseudo-Sculling Errors 


Table 5.2-2 summarizes the sensor-level and system-level 
errors induced by the helicopter environment. This table has been pre- 
pared in the same format as the error summary for a system based on 
the Norden 1139 single-degree-of-freedom gyro (Table 4.2-3, Ref. 4) 
in order to facilitate a better understanding of the relative capabilities 
of the two types of gyros. The rebalance loop asymmetry error for the 
SDG-1 gyro is about half the error obtained in Ref. 4 for the Norden 113 
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TABLE 5,2-2 

VIBRATION INDUCED ERRORS: SDG-1 GYRO 


Vehicle: 

CH-46C Twin Rotor Helicopter 

Gyro: 

Teledyne SDG-1 Gyro 

Accelerometer: 

Perfect Accelerometer Assumed 

Principal Sources of Vibration -Induced Error 

Gyro Drift Rates <deg/hr) 

Axes 

x y z 


• Anisoinertia 

-0.001 

0.009 

-0.003' 

. • Cross. Coupling 

-0.011 

0.091 

0:081 

• Angular Acceleration Sensitivity 

— 

0.001 

0.001 

• Rebalance Loop Asymmetry 

0.065 

0.061 

0.035 

System Drift Rates (deg/hr) 




* Angular Acceleration Sensitivity 

0.075 

0.003 

0.201 

• Undetected Coning Motion 

0.042 

-0.115 

0.310 

System Acceleration Errors (yg) 




• Angular Acceleration Sensitivity 

0.28 

0.40 

-0.11 

• Undetected Sculling Motion 

-15.5 

2.8 

8.1 


gyro. This is a direct consequence of the fact that the asymmetry 
parameter used for the SDG-1 was half that used for the Norden 1139. 

The most significant system-level error mechanism is the angular accel- 
eration sensitivity exhibited by the SDG-1. This produces pseudo-coning 
errors at the system level that are of the same order of magnitude as 
those produced by a similar effect in the Norden 1139. 

The undetected coning motion error is somewhat larger for the 
SDG-1 than for the Norden 1139: Recall that undetected coning motion is 
coning motion that occurs at a frequency above the detection bandwidth of 
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the gyro. The SDG-1 rebalance loop bandwidth used in computing the - 
errors listed in Table 5.2-2 was 25 Hz while the bandwidth for the 
Norden 1139 was 70 Hz. This explains the difference in the undetected 
caning motion errors. 

The SEPP program was also exercised'for a* system contain- 
ing SDG-1 gyros in which the bandwidth of the gyros had'been increased 
to about 50 Hz, through appropriate changes in rebalance; loop gains and 
compensation. For purposes of comparison the major^errors* generated 
in that system are shown in Table 5.2-3. Tables 5.2-2 and 5.2-3 reveal 
that the gyro loop with a wider bandwidth has significantly smaller errors 
due to cross coupling and undetected coning and sculling motion. The 
cross coupling errors are reduced because the gyro has a stiffer rebalance 
loop and the undetected coning and sculling are less because the gyro cor- 
rectly measures more of the motion spectrum. Comparing the two tables 
also shows that rebalance loop asymmetry errors go up, also because 
more of the motion spectrum is being passed by the gyro. Pseudo-coning 
errors caused by angular acceleration sensitivity are lower in Table 5.2-3, 
evidently because of the tighter rebalance loop. Of course, certain errors 
not related to baixi width, such as anisoinertia, are unchanged. The 
total of the major deterministic errors (i.e. , excluding those with ran- 
dom error coefficients such as rebalance loop asymmetry and rotor 
bearing compliance) is somewhat smaller for the gyro with the 50 Hz 
bandwidth, mainly because of the reduction in z-axis undetected coning. 

It is not clear that the larger bandwidth is generally superior for all 
motion environments. 
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TWt.NALTTK SCUMCES CORPORATION 


TABLE 5.2-3 

VIBRATION INDUCED ERRORS: 
SDG-1 GYRO WITH WIDE BANDWIDTH 


Vehicle: CH-46C Twin Rotor Helicopter 

Gyro: Teiedyne SDG-1 Gyro (Modified) 

Accelerometer: Perfect Accelerometer Assumed 

Principal Sources of Vibration -Induced Error 

Gyro Drift Rates (deg/hr) 

Axes 




51 

z 

• 

Anisoinertia -0.001 

0.009 

-0.003 

• 

Cross Coupling -0.001 

— 

-0.001 

• 

Angular Acceleration Sensitivity 

— 

— 

6 

Rebalance Loop Asymmetry 0.074 

0.069 

0.013 

System Drift Rates (dec/hr) 



• 

Angular Acceleration Sensitivity -0.067 

0.002 

0.168 

e 

Undetected Con mg Motion — 

— 

— 

System Acceleration Errors (wr) 



• 

Angular Acceleration Sensitivity 0.33 

0.24 

-0.11 

0 

Undetected Sculling Motion -2.3 

-3.1 

2.4 


5.3 SUMMARY 

The error propagation characteristics of a navigation system 
comprised of two Teledyne SDG-1 gyros and three perfect accelerometers 
has been studied using the TASC System Error Prediction Program. The 
dynamic motion environment of the CH-46C helicopter was used as the 
excitation function. The major error mechanisms and the magnitudes of 
the errors produced by these mechanisms are summarized in Tables 
5.2-2 and 5.2-3 for two different sets of rebalance loop gains and com- 
pensations. 


6. SUMMARY AND CONCLUSIONS 

This report has described the operation of a tuned, flexure - 
mounted gyro as a strapdown sensor* Transfer functions describing the 
gyro dynamic performance are derived and discussed for the case of a 
linear rebalance loop. 



A dynamic error model is derived for the gyro, describing the 
instrument's sensitivity to 


• Angular, motion 

• sipecihc force" 

• Rotor speed variations 



" h order to obtain an -understanding of the performance capa- 

bility of a strapdown navigation system utilizing a pair of tuned, flexure - 
mounted gyros, a sample case was investigated using the TASC System 
Error Propagation Program. The Teiedyne SDG-1 gyro was chosen as 
the system gyro and the dynamic motion environment was assumed to be 
tha^tf^a CH-46 helicopter. Leading contributors to individual gyro drift 
were found to be : 


• Cross coupling 

• Rebalance loop scale factor asymmetry 


The dominant contributors to system-level errors such as coordinate 
frame drift were : 


5-10 


6-1 



IE ANALYTIC SCIENCES CORPORATION. 



• Sensitivity to angular accelerations 

• Bandwidth limitations on the gyro 


In general a system containing the SDG-1 gyro was found to have dynamic 
errors comparable in magnitude to a system (investigated in Ref. 4) which 
employs the Notrden 1139 gyro. 


A 'second set of calculations were made for afcSDG-1 gyro with . 

■ j j ■ •• 

an increased bandwidth. The system errors calculated were smaller , ^but 

■ • * * ' ‘ . ‘ Jr 1 ■ * y * 

nnl' cionifirantlfT CG 
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APPENDIX A 

ANGULAR MOTION TORQUES 

In this appendix we compute the rotor's reaction torque due to 
gyro case angular velocity. The reaction torque vector L is given by 

' L = [H], (A-l) 

~ ~ 1 . 

where H is the rotor angular momentum vector and [ ]\ indicates that the 
derivative is taken with respect to an inertial coordinate frame. Apply- 
ing the Law of Coriolis to this relationship yields the torque in terms of 
derivatives relative to the rotor -fixed coordinate frame 

L = [H] r + to x H (A— 2 ) 

where is the angular velocity of the rotor relative to an inertial frame 
and [ ] r indicates that the derivative is taken with respect to a coordinate 
frame related to the rotor principal axes. In particular the coordinate 
frame we will use has two axes in. the plane normal to the rotor axis of 
symmetry and one axis along the rotor spin (symmetry) axis. The first 
two axes are not fixed in the rotor, however. They rotate with respect 
to the rotor at a speed equal to and opposite in direction from rotor spin. 
The w vector is the sum of the gyro case angular velocity vector G and 
the gimbal angle derivative vector <5. 

u = O + 0 (A-3) 


A-l 


• 




THE ANALYTIC SCIENCES CORPORATION 


In gyro case coordinates, 



°x + V 
V*y 

°z - 


(A-4) 


The coordinate frame most natural for Eq. (A-2) is the rotor referenced 
coordinate frame* since the axes are along principal inertia axes of the 
rotor. 


i 


To convert from case to rotor coordinates, we apply the co- 
ordinate transformation, which for small rotor deflection angles 0 is 
given by Eq. (2.1-1): 



Applying this to Eq. (A-4) and neglecting second-order products of rotor 
deflection angles and derivatives, we obtain 


’n +e *8fi 

X x y z 

( w) = 0+0+80 

r y y x z 

n o - e o 

z y x x y 

* 

As explained in Chapter 2, the third (C) axis of the rotor referenced 
frame is parallel to the rotor axis of symmetry although this frame is 
non-rotating relative to the gyro case about this axis. For convenience, 
the rotor referenced frame is referred to as the rotor frame. 
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The angular momentum vector H is 

H = j(&>+ w s ) (A -7) 

where for rotor principal axes, the inertia tensor J is 

AO O' 

0 A 0 (A-8) 

0 0 C. 

and the rotor spin angular velocity oj g is, in rotor axes, 

0 
0 

w s 

Combining the four previous equations, the angular momentum, resolved 
in the rotor coordinate frame, is 

Afn + 6 -en) 

V x x y z/ 

(2)r = A (°y +6 y + 6 x° Z ) (A-10) 

c fo 7 + e„n - e o + «) 

V z y x x y s/ 

By inserting Eqs. (A-6) and (A-10) into Eq. (A-2), we obtain in 
the rotor 5, »», C coordinates shown in Fig. 2.1-1, 

L i = A (V S x) +H (Vy +(C - A) 0 y°z 

+ 9 x[( C - A )(°z-°y) + HO z ] 

+ e y [(c-A)n y o x - A n z ] (A-ii) 
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Then, using Eqs. (A-14) and (A-5), we get 


j L = L, + 6 L, 

(A-17) 

. i X 1 y c 


1 L = L - 8 L 

(A-18) 

i i y v x c 



The z component of the torque vector is of no interest since it is; about 
the rotor shaft axis and not a 0 x or 8^ rotor deflection pick-off axis. 
Using Eqs. (A-11) through (A-13) in Eqs. (A-17) and (A-18), we have 

’ L x = A ( rt x + S x) +H (°y + 6 y) + < C - A )°y b z 
+ 8 x [ (C ' A) (°z ‘ °y ) +Hn z ] 

+ e y [(c-a) n y n x + (c-a) A J (a-w) 

and 

L y = A ( A y + § y)- H ( n x + 0- (C - A )°x n z 

• ' e y [(C' A )( n z'°x)' HO z] 

+ e x [(c-A)nn x -(c-A)nJ (A-20) 

We may identify the terms HO in L and -HO in L as the "ideal" tor- 

y x , x y 

que components, i. e. , the components which in an ideal gyro would be 
attributed to the case angular motion. The other terms contribute to 
gyro errors as is shown in Chapters 3 and 4. 
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APPENDIX B 

SPECIF k: FORCE TORQUES 

Two types of error torques are considered in this Appendix to 
arise from the specific force experienced by the gyro case. They are the 
mass unbalance torques , given by 

ij = m(6xf) ’ (B-l) 

where 

m is the rotor mass 

is the displacement of the center of mass from 
the center of support 

i is specific force 

and the compliance torque, given by 

£ 2 = m 2 (f x Kf) (B-2) 

where K is the rotor compliance matrix. 

For the mass unbalance torque, we will take 6 to be, in case 
coordinates, 

O' 

0 (B-3) 

6 

z 
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since the rotor is spinning at a virtually constant angular velocity and 
the z component alone has a non -zero average value. Likewise in case 
coordinates, 


Using Eqs. (B-3) and (B-4) in Eq. (B-l) yields for the x and y components 


For the compliance torque, the most general form of the com- 
pliance matrix K is 


|K K K 1 
xx xy xz I 


K K yx K yy K yz 
K zx K zy K zz 


Using this and Eq. (B-4) in Eq. (B-2) yields for the x and y components 


•2, k«y, * v? * - va - v*v - vJ 

= m 2 

. * 2 y . K xx f z f x * K xy f z f y + K xz*z “ K zx*x " K zy f x f y " K zz Vz 
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